Polarity patterns of stress fibers 
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Stress fibers are contractile actomyosin bundles commonly observed in the cytoskeleton of meta- 
zoan cells. The spatial profile of the polarity of actin filaments inside contractile actomyosin bundles 
is either monotonic (graded) or periodic (alternating). In the framework of linear irreversible ther- 
modynamics, we write the constitutive equations for a polar, active, elastic one-dimensional medium. 
An analysis of the resulting equations for the dynamics of polarity shows that the transition from 
graded to alternating polarity patterns is a nonequilibrium Lifshitz point. Active contractility is a 
necessary condition for the emergence of sarcomeric, alternating polarity patterns. 

PACS numbers: 47.54.Bd, 47.54.Fj, 87.16. Ka, 87.16. Ln 



Stress fibers are bundles of actin filaments made con- 
tractile by interaction with myosin minifilaments and 
elastic by crosslinking with a-actinin and other proteins. 
They are formed in the cytoskeleton of nonmuscle animal 
cells that need either to exert or to resist a mechanical 
force (see pQ for a recent review). Their biological func- 
tion is exemplified in vivo by myofibroblasts that use the 
contractility of stress fibers to remodel the extracellular 
matrix during wound healing 2J , or by vascular endothe- 
lial cells in blood vessels that align stress fibers parallel 
to the direction of flow to resist shear forces [3] . 

Actin filaments are polar structures, with a plus (or 
barbed) end and a minus (or pointed) end. The polarity 
of actin filaments inside a bundle is determined by dec- 
oration of actin with myosin sub fragment- 1 and electron 
microscopy imaging. Noncontractile (or passive) actin 
bundles often present uniform polarity, with barbed ends 
predominantly facing the same direction along the bun- 
dle, as seen in filopodia [3] or in the retraction fibers of 
postmitotic cells [S]. The polarity profile of stress fibers 
is alternating, with a majority of barbed ends succes- 
sively pointing in opposite directions O Ej • This peri- 
odic structure, with a wavelength of the order of 1 fim, 
is reminiscent of the sarcomeric organization of muscle 
cells. The actomyosin filament bundles of motile fibrob- 
lasts, however, exhibit a graded polarity profile [8]: 
The proportion of barbed ends in the bundle facing one 
direction varies gradually and monotonically between 
and 100 % from one end to the other. The direction of 
the bundle sets the orientation of locomotion The 
formation of graded polarity bundles is an early event in 
the polarization of fibroblasts prior to migration. The 
spatial dependence of polarity in contractile actomyosin 
bundles contrasts with the uniform polarity observed in 
structures such as cilia, filopodia, and retraction fibers, 
in which motors do not play an organizing role. 

In this Letter, we aim at identifying the physical mech- 
anisms responsible for the existence of variable polarity 
patterns [10] in elastic [HI [12] actin bundles. We show 
that active contractility, the chcmomcchanical transduc- 
tion from ATP hydrolysis to motion by myosin motors 



"pulling" on actin, can generate periodic, sarcomerelike 
polarity patterns provided that it exceeds a well defined 
threshold. We focus our attention on ventral stress fibers 
PQ that span the cell ventral surface and transmit force 
to the underlying substrate through end-point focal ad- 
hesions. We consider existing actomyosin bundles, with 
constant total material content, and ignore the intriguing 
questions raised by their assembly |13j . 

Recent theoretical descriptions of stress fibers focus on 
their mechanical properties, e.g., viscoelastic recoil after 
photoablation [T4] . The question of polarity sorting was 
addressed in early models of contractile bundles of actin 
filaments [TS] that consider a liquid state in the (dilute) 
limit of low concentration. In a recent extension, density 
patterns for viscoelastic actomyosin bundles are studied 
in the same limit [16] . The effect of coupling the density 
and polarity fields in the liquid state is investigated in 
detail in Ref. [17]. Conditions for the spontaneous oscil- 
lations of muscle sarcomeres are determined within the 
context of active elasticity in Ref. [TB]. We propose to 
formulate the constitutive equations for an active, polar, 
elastic one-dimensional medium, following the prescrip- 
tions of linear irreversible thermodynamics [THl HO] i valid 
in principle for arbitrarily large densities. This work is 
directly inspired by Ref. [5T], where constitutive equa- 
tions were proposed for an active polar viscoelastic liquid 
(see [22] for a review). Active polar liquid crystals share 
the same invariance properties j23M25]. We treat here 
the case of a one-dimensional viscoelastic solid (Kelvin 



actin filament (+1) 


• a-actinin 


•— • myosin filament 


|r* focal adhesion 




FIG. 1. (Color online) Schematic representation of an ac- 
tomyosin bundle with graded (A) and periodic (B) filament 
polarity profile. 



2 



model), with an emphasis on how polar degrees of free- 
dom couple to elastic strain. Within this formalism, we 
restrict our attention to lowest-order terms in the free 
energy and in the constitutive equations. 

In the spirit of Ref. [3T], we model a stress fiber as a 
one-component material. The description is mesoscopic, 
valid over length scales large compared to the typical 
mesh size of the polymer network (£ ~ 10 nm). The 
geometry is one-dimensional: The stress fiber of (con- 
stant) length L occupies the interval x <E [-L/2, L/2]. 
Our hydrodynamic description is built on three conserved 
fields: the mass, momentum and energy densities p(x,t), 
g(x,t), and u(x,t), respectively; and one slow mode, the 
polarity p — p(x,t) e x , defined as the coarse-grained av- 
erage of the polarity of actin filaments present at time t 
in the vicinity of x. By convention, p = +1 (respectively, 
p = — 1) corresponds to a bundle where all (respectively, 
no) barbed ends point towards positive x (see Fig. [I]). 
The relevant boundary conditions for the polarity field 
are "antisymmetric" : 



p(-i/2) = -l, p(i/2) = +l 



(1) 



since the barbed ends of actin filaments face the end- 
point focal adhesions [H [7] . 

The stress fiber is a viscoelastic solid [12] , with velocity 
and strain fields v(x,t) and e(x,t). We focus on time 
scales for polarity ordering that are long compared to the 
typical turnover time of actin and other dynamic protein 
components: t 3> tto — 10 2 s [15] . The density field p is 
a fast variable, slaved to the strain field: Sp/p = — e. The 
free energy density, expanded to quadratic order about 
p = e = 0, reads [2"5] : 

/=^ 2 + f (d xP ) 2 + ~e 2 +wed xP , (2) 

where all terms are invariant under inversion of space 
and polarity. The positive coefficients a, K, and G re- 
spectively denote the susceptibility, the energy cost of in- 
homogeneities of the polarity field, and the elastic mod- 
ulus. Symmetries of the problem allow for a coupling 
term between strain and thepolarity gradient, with cou- 
pling parameter w of unknown sign. From Eq. j2j), we 
obtain the molecular field h, conjugate to the polarity: 
h = — = — a p + K d x p + w d x e and the elastic 

stress cr el = = G e + w d x p. Thermodynamic stability 
imposes the constraint w 2 < K G: The uniform, strain- 
free, zero polarity steady-state is a stable solution when 
boundary conditions permit. 

We proceed to derive the dissipation rate R: 



R=(a + P- <j c1 ) d x v+ph + r Ap, 



(3) 



where P and a denote the pressure and (total) stress 
fields, respectively. Equation ^ contains a mechanical 
term, a polar term {p is the total derivative of p), and 



an active term where the chemical force Ap is conjugate 
to the ATP consumption rate r [21 . We next express 
the relevant generalized thermodynamic fluxes as linear 
functions of the generalized thermodynamic forces and 
obtain the stress fiber's constitutive equations: 



P 



P 



r h - a pd x p Ap (4) 
V d x v + (-( + p d xP ) Ap (5) 



The positive kinetic coefficient T is analogous to an in- 
verse viscosity [37], and rj is the viscosity of the fiber. A 
polar, nonrelaxational term with a coupling coefficient a 
of unknown sign is included in Eq. Q, in agreement with 
studies of active liquid crystals [231425] . The last term in 
Eq. ([5|, a nondiagonal chemomechanical coupling, gives 
the active contribution to the total stress. The coeffi- 
cient C is negative for contractile actomyosin 21 . The 
invariance properties of polar media allow for an active 
coupling between stress and polarity gradient, with a co- 
efficient /3, as first introduced in Ref. [55] • For simplic- 
ity, we ignore the punctate spatial pattern of myosins in 
stress fibers Q3[Z] : A/i is a constant in the following. 

We further assume that the time scale for polarity or- 
dering is much longer than that for strain relaxation, of 
the order of 1 s as measured in laser ablation experi- 
ments [12] . In this long time limit the velocity field 
v and its gradient are negligible. The pressure field 
is uniform and equal to the pressure in the surround- 
ing cytosol P = Pq. Newton's law, d x a — 0, yields: 
G d x e + (w + (3 Ap) d x p = 0, which eliminates the strain 
field e in Eq. Q. This yields the evolution equation for 
the polarity field. 

In the passive case Ap = 0, the total derivative p re- 
duces to d t p when v = 0. We obtain a damped diffusion 
equation for the polarity field: 



d tP = -aT (p - I 2 d 2 xP ) 



(6) 



where the length scale I 2 = (K/a) [l - w 2 /(KG)] > 
has been introduced. Using antisymmetric boundary 
conditions, we find a linearly stable, stationary, mono- 
tonic solution of Eq. p(x) = sinh(x/l) / smh(L/2l) 
that describes the graded polarity profile of actin bundles 
in the absence of myosin-induced contractility. 

In the active case Ap > 0, Eq. Q yields a damped 
Burgers equation: 



dtp + aAp pd x p = —aTp + D d 2 x p, 



with a diffusion constant D = TK 



1 



1 



KG "T" w J 

For weak activity, i.e. when (3Ap/w < KG/w 2 — 1, D 
is positive, and monotonic stationary solutions qualita- 
tively similar to the above analytical expression are ob- 
tained numerically (see Fig. [2]). These linearly stable 
solutions describe the graded polarity patterns observed 
in contractile stress fibers. 
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FIG. 2. (Color online) Stationary polarity profiles obtained 
from numerical simulations. We set the dimensionless system 
size to L = 128 and use the boundary conditions Top: 
When ii = and D = 500, the graded profile obtained for 
the damped Burgers equation (black line) agrees qualitatively 
with experimental data, taken from Cramer et al. [7], where 
the first and last data points are set at the boundary (circles). 
Mixed polarity is obtained far from the tips when D = 10 
(red line). Bottom: Typical periodic profile for the damped 
KS equation with v — 1.0 and D — —2.01 (black line). Mixed 
polarity is obtained far from the tips with v — 1.0 and D = 
-1.90 (red line). 



For strong activity, i.e. when /3A/i/w > KG/w 2 — 1, 
the diffusion coefficient changes sign to become negative. 
In the particular case a = 0, we obtain periodic station- 
ary solutions: p(x) — sin(qx) / sm(qL/2), with a wave 
number q 2 = —aT/D. Since D < 0, stationary solu- 
tions of Eq. ([7]) are linearly unstable. We are led to in- 
clude a stabilizing, higher-order polarity gradient term 
vj2 [d 2 p) in the free energy density, where v is a pos- 
itive higher-order diffusion constant. This is reminiscent 
of a Lifshitz point |20[ 12 7j , where an equilibrium phase 
transition between a uniform and a modulated phase oc- 
curs as a diffusion coefficient changes sign. The field h is 
supplemented with the term — v d x p and Eq. Q becomes 
a damped Kuramoto-Sivashinsky (KS) equation: 

d t p + a Afi pd x p =-aTp + D d 2 x p - vT d^p, (8) 



For large enough damping (here aT), this equation pos- 
sesses stable stationary solutions periodic in space |28j . 
The transformation t = aTt, x = aTx/{\a\ A//), makes 
Eq. (Isj) dimensionless, with only two parameters D = 
( aA M )2 and ^ = a(^A7l) 4 - Alternating polarity patterns 
are obtained when the uniform solution p{x) = is lin- 
early unstable, i.e. for D < D c = —2*/i> (see Fig. 
Their wavelength A is obtained where the growth rate of 
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FIG. 3. Phase diagram of an active fiber in the (D, ii) plane, 

The polarity profile is 
> +00 when A/i — > 0). 
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the linearized damped KS equation is maximal: 
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(9) 



We predict that the inverse square wavelength A~ 2 of al- 
ternating polarity patterns is a linear function of A/i, 
as may be checked experimentally by up- or down- 
regulating the contractility of fibers. The wavelength 
increases with G and is independent of cv, and also of 
L when L ^> A. We checked numerically that the pe- 
riodicity of patterns is robust to sufficiently weak spa- 
tial modulations of Afi mimicking the punctate profile 
of myosins, as well as to the addition of a noise term of 
small amplitude. 

The case of "mixed" polarity (p = 0) is also ac- 
counted for [TU]. When D > 0, the uniform solution 
p(x) = e{x) = is stable but irrelevant because of in- 
compatible boundary conditions. Mixed polarity may, 
however, be observed locally, far from the tips, when the 
penetration length 8 cx \/~D is small compared to the to- 
tal length L (see Fig. |2j. When D < 0, mixed polarity is 
also obtained far from the tips in the range D c < D < 0, 
with alternating polarity near the boundaries. Simula- 
tions suggest that the penetration length is of the order 
of S oc (D - D c )~ 1/2 . To the best of our knowledge, 
this type of pattern has not been observed in experi- 
ments. These results are summarized in a phase diagram 
(Fig.§. 

In the above derivation, we implicitly assumed that 
polarity is a nonconserved field. By disregarding the 
possibility of spontaneous polarity reversal of a filament, 
which we deem extremely unlikely in a one-dimensional 
cross-linked bundle, this implies that relevant time scales 
are longer than the characteristic time tq for nucle- 
ation/annihilation and insertion/removal of actin fila- 
ments. Equations Q and ([8]) are valid only when t > tq. 

We now derive the constitutive equations for a con- 
served polarity field, valid for shorter time scales t < tq 
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(see [T7] for an alternative approach). The conservation 
equation reads p = —d x j p , and j p , the polarity current, 
is the flux conjugate to the force d x h. The constitutive 
equations read: 

f = Td x h+ (-C + pd x p+a p 2 ) A/i - X d x v 
cr = -P + a cl + V d x v + (-( + pd x p)Ati + X d x h 

where £, f3, and x are additional cross-coupling coeffi- 
cients of unknown sign and the nonlinear term ap 2 Afi 
is included in analogy with the nonconserved case. An 
Onsager relation sets the reactive nondiagonal coupling 
between stress and the gradient of the molecular held. 

For a passive bundle (A/i = 0), we obtain the dif- 
ferential equation obeyed by stationary polarity profiles: 
d x (p — I 2 d 2 p) = 0, in the limit where v = d x v — 0. Lin- 
ear stability of the uniform, zero polarity steady-state im- 
poses the condition wx > 0. For antisymmetric bound- 
ary conditions [Eq. Q], we find as a linearly stable sta- 
tionary solution the linear profile p(x) — 2x/L. 

For an active bundle (A/i > 0), we obtain a set of two 
coupled equations for p and h: 



d t p + aA/j p d x p = -T d 2 h - /?A/i d 2 x p, 



(10) 



wx 
G 



dl)h 



D 



ap + ^ d x p - v d x p, (11) 



where the higher-order term v/2 {d 2 p) was included in 
the free energy density. An exhaustive study of the phase 
diagram of the system ( TofTT I is beyond the scope of this 
Letter. We naturally expect the graded polarity patterns 
obtained in the passive case to be preserved for weak 
enough activity. We shall restrict our analysis to the limit 
wx » G, where the molecular field may be eliminated. 
In this limit, we obtain a damped KS equation: 

d t p + aA/i pd x p = -e c p + D c dip - v c d x p (12) 



with a damping coefficient e c 



GTa/wx, a diffu- 
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sion constant D c given by: 

A/i \~x~^P)' anc ^ a higher-order diffusion constant 
v c = GTu/wx- In this regime, beyond the nonequilib- 
rium Lifshitz point (D c < 0), and for large enough damp- 
ing e c , we shall again find alternating polarity patterns. 
The transition may be driven either by the active stress 
(coefficient j3) or by the active current (coefficient /3). 

A key ingredient of our model of an active, polar, one- 
dimensional elastomer is the coupling between the elastic 
strain and polarity gradient, allowed by symmetry in the 
free energy expansion [Eq. |[2j)]. We show that an al- 
ternating polarity profile emerges in an active medium 
with a uniform spatial distribution of myosins: Cou- 
pling with the myosin density field is expected to lead to 
myosin patterns with the same wavelength. While graded 
polarity patterns are allowed for both passive and active 
bundles, active contractility is a necessary condition for 



the self-organization of stress fibers into sarcomerelike, 
alternating polarity patterns. The alternating polarity 
pattern of an active actomyosin bundle may therefore 
become graded upon treatment with drugs that inhibit 
contractility. Whereas alternating profiles occur far from 
the boundaries beyond the instability threshold, graded 
profiles are made possible by nonzero values of the polar- 
ity field at the tips. These conclusions hold irrespective 
of the conserved or nonconserved nature of the polarity 
field. 

Linear irreversible thermodynamics is a powerful fra- 
mework within which coupling terms are deemed relevant 
according to their invariance properties. The numerical 
values of coupling constants are beyond the scope of the 
theory and may well depend here on the cell type. The 
only measurement we are aware of is that of the elas- 
tic modulus G ~ 10 5 — 10 6 Pa Giving a micro- 
scopic interpretation to hydrodynamic coefficients such 
as w, which couples elastic and polar degrees of freedom, 
will require models that relate behavior at hydrodynamic 
scales with the microscopic interaction of actin filaments 
with active and passive cross-linkers. One possible in- 
terpretation of w is that a preference of cross-linkers to 
parallel or antiparallel pairs of actin filaments may lead 
to a different value of the stress at constant strain. 



t 
t 

[1 

[2 
[3 

K 

[5 

[6 
[7 
[8 
[9 
[10 

[H 
[12 

[13 

[14 
[15 



[16 
[17 
[18 

[19 
[20] 



Present address: Fukui Institute for Fundamental Chem- 
istry, Kyoto University, Kyoto 606-8103, Japan 
E.S.P.C.I., 10 rue Vauqu elin, 75005 Paris, France 
philippe.marcq@curie.fr 

S. Pellegrin and H. Mellor, J. Cell Sci. 120, 3491 (2007). 
J. Tomasek et al, Nat. Rev. Cell Mol. Biol. 3, 349 (2002). 
D. Drenckhahn and J. Wagner, J. Cell Biol. 102, 1738 
(1986). 

A. Verkhovsky et al, J. Cell Sci. 110, 1693 (1997). 

L. P. Cramer and T. J. Mitchison, J. Cell Biol. 131, 179 

(1995). 

J. Sanger and J. Sanger, J. Cell Biol. 86, 568 (1980). 

L. P. Cramer et al, J. Cell Biol. 136, 1287 (1997). 

T. Mseka et al, Cell Motil. Cytoskeleton 66, 743 (2009). 

T. Mseka et al, J. Cell Sci. 120, 4332 (2007). 

L. P. Cramer, Biochem. Soc. Symp. 65, 173 (1999). 

S. Deguchi et al, J. Biomech. 39, 2603 (2006). 

S. Kumar et al, Biophys. J. 90, 3762 (2006). 

P. Hotulainen and P. Lappalainen, J. Cell Biol. 173, 383 

(2006). 

J. Colombelli et al, J. Cell Sci. 122, 1665 (2009). 
H. Nakazawa and K. Sekimoto, J. Phys. Soc. Japan 65, 
2404 (1996); K. Kruse and F. Jiilicher, Phys. Rev. E 67, 
051913 (2003). 

R. Peter et al, New J. Phys. 10, 035002 (2008). 
K. Kruse et al, Europhys. Lett. 64, 716 (2003). 
S. Giinther and K. Kruse, New J. Phys. 9, 417 (2007). 
S. R. de Groot and P. Mazur, Non-Equilibrium Thermo- 
dynamics (Dover Publications, 1985). 
P. Chaikin and T. Lubensky, Principles of Condensed 
Matter Physics (Cambridge University Press, 2000). 



5 



[21] K. Kruse et ai, Eur. Phys. J. E 16, 5 (2005). 

[22] F. Jiilicher et ai, Phys. Rep. 449, 3 (2007). 

[23] A. Ahmadi et ai, Phys. Rev. E 74, 061913 (2006). 

[24] S. Muhuri et ai, Europhys. Lett. 78, 48002 (2007). 

[25] L. Giomi et ai, Phys. Rev. Lett. 101, 198101 (2008). 

[26] Since the boundary conditions on p are fixed, we may 



ignore the polar term d x p in the free energy density. 

[27] P.-G. de Gennes and J. Prost, The Physics of Liquid 
Crystals (Oxford University Press, 1995). 

[28] P. Manneville, Dissipative Structure and Weak Turbu- 
lence (Academic Press, 1990). 



